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The Ludvigsen-Vickers and two recently suggested quasi-local spin-angular momentum expressions, based 
on holomorphic and anti-holomorphic spinor fields, are calculated for small spheres of radius r about a 
point o. It is shown that, apart from the sign in the case of anti-holomorphic spinors in non-vacuum, the 
leading terms of all these expressions coincide. In non-vacuum spacetimes this common leading term is of 
order r 4 , and it is the product of the contraction of the energy-momentum tensor and an average of the 
approximate boost-rotation Killing vector that vanishes at o and of the 3-volume of the ball of radius r. In 
vacuum spacetimes the leading term is of order r 6 , and the factor of proportionality is the contraction of the 
Bel-Robinson tensor and an other average of the same approximate boost-rotation Killing vector. 

1. Introduction 

In the classical theory of matter fields the energy-momentum and angular momentum density are described 
by the symmetric energy-momentum tensor T ab . For Killing vectors K a the current T ab K}, is conserved, and 
it is interpreted as a component of the energy-momentum or angular momentum current density, depending 
on the nature of the Killing field. E.g. the translation- and rotation Killing vectors (or rather 1-forms) 
of the Minkowski spacetime are given by :— V e x a and by := x a K^ — x h Kf, where x a are the 
standard Descartes coordinates, and the corresponding currents are the four conserved energy-momentum 
and six conserved total (i.e. orbital- and spin-) angular momentum currents, respectively. By converting the 
name indices of the Killing 1-forms to spinor name indices the translations form a constant Hermitian matrix 
valued 1-form K AA , and the rotations can be decomposed as K AA BB = e A B K AB + e AB K A B , the sum 
of the anti-self-dual and self-dual rotations, where e.g. the anti-self-dual rotations are given explicitly by 
K AB :=x( A B ,K? )B ' . If S is any smooth spacelike hypersurface with fixed smooth 2-boundary $ := <9£, t a 
is its future directed unit normal and d£ is the induced volume element, then P AB := J s T ab t a K AB dS 
and J AB :— J E T ab t a K AB AT, depend only on the boundary $ (and are independent of the actual £ defining 
the homology between $ and zero) , and hence may be interpreted as the quasi- local energy-momentum and 
angular momentum (with respect to the origin o of the Descartes coordinate system) of the matter fields 
associated with $, respectively. Let r be a positive number and £ := {(t,x,y, z)\t = r > \J x 2 + y 2 + z 2 }, 
i.e. the piece of the spacelike hyperplane t = r that is bounded by its intersection with the future null cone 
Af„ of the origin o, and let $ r := 9S. Then the quasi-local energy- momentum and angular momentum can 
be computed for the sphere $ r of radius r by taking the flux integrals of the conserved currents on S. But 
since P AB and J AB depend only on $ r and the origin o is a regular point, they equal to the flux integral 
of the conserved currents on the null cone Af Q between the vertex o and $ r . (For the technical details, e.g. 
the volume element on the null N a , see Section 3 below.) For sufficiently small r the quasi- local energy- 
momentum is r-Trr 3 {T ab t a K^ B ')\ , the product of the 3-volume of S and the value of T ab t a K AB ' at the 
origin. To compute the angular momentum too, it seems useful to introduce spinors. Let := {0 A ,I A } 
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be the constant normalized spin frame field associated with the Descartes coordinates, £ A := {—Ia, Oa} its 
dual basis, and e\ := {o A , i A }, := {— la, oa} the pair of standard dual spin frames on the null cone Af - 
If £ :— expi(/>cot |, the complex stereographic coordinates (and hence r 7 (,( form a coordinate system on 
Mo - {o}), then o A (r, (, () = i^2(l + ((O a + I A ) =: X A £ A , and hence X A = £%o A . On the null cone 

the Hcrmitian matrix :r AB of the Descartes coordinates becomes r-times the dyadic product of the 'spinor 
coordinates': x AB = rX A X B . Thus on the null cone the anti-self-dual rotation Killing field becomes 

r/ AB Ac(A c B)_ -i \ 

K e = ro £\ £ E 'o E >, (1-1) 
and, with accuracy r 4 , the quasi-local anti-self-dual angular momentum for small spheres of radius r is 

jAB = y^ b £^£^)\ ^OAO A ^-o B ,A% = ^T A A'BB't AA 't B ' E e BF 4 A £P\ . (1.2) 

Here $ is the unit sphere and d$ :— — 2i(l + ££)~ 2 d(d£, the 2-surface element on $. Therefore the quasi-local 
angular momentum for small sphere of radius r is -|r 3 times the contraction of the energy- momentum tensor 
at o and an average of the anti-self-dual rotation Killing vector that vanishes at o. 

As a consequence of the diffeomorphism invariance of the geometric theories of gravity every vector 
field generates a conserved current (and its flux a conserved quantity), which can always be derived from 
a globally defined supcrpotcntial 2-form, independently of the homological structure of the spacetime [1]. 
In the case of energy-momentum in Einstein's theory this superpotcntial appears to be the Nester-Witten 
2-form u(\,fi) a b, associated with any pair of spinor fields Xa and \xa [2,3], because of the following reasons. 
First, this 2-form defines a 2-form on the spin frame bundle, and hence on the bundle of orthonormal 
frames over the spacetime manifold too, which 2-form extends uniquely to the bundle of linear frames L(M). 
The superpotentials of the various classical energy-momentum pseudotensors (e.g. the Einstein, Bergmann, 
Landau-Lifshitz and the tetrad-M0ller pseudotensors) are just the pull backs of various forms of this 2-form 
along various local cross sections of L(M) [4,5]. In fact, the exterior derivative of the Nester-Witten form, 
known as Sparling's equation, looks like the Noether identity: V[ a u(A, p)bc\ — T(A, fi)abc — ^X£,flD'G de ^e ea i,c 7 
where G a b is the Einstein tensor and I\A, p,) a bc, the so-called Sparling 3-form, is a quadratic expression of 
the derivatives of the spinor fields [6,7]. Second, both the ADM and Bondi-Sachs four-momenta can be 
written as the 2-surface integral of the Nester-Witten 2-form at spacelike and null infinity, where the spinor 
fields are chosen to be the constant or the asymptotic spinors there, respectively. These spinors may also 
be interpreted as the spinor constituents of the asymptotic translations at infinity. The simplest proof of 
the positivity of the ADM and Bondi-Sachs masses is probably that based on the use of the Nester-Witten 
2-form [2,3,8-13]. Finally, the integral of the Nester-Witten 2-form for a spacelike topological 2-sphere $ in 
the spacetime can be used to define energy-momentum, a manifest Lorentz covariant quantity (and not only 
energy, or energy and linear momentum separately), at the quasi- local level too. The only question is how 
to choose the two spinor fields Xa and \ia- Ludvigsen and Vickers [14] suggested a rule of transportation of 
the asymptotic spinors from infinity back to $ if the spacetime is asymptoically flat at future null infinity 
and $ can be connected with the future null infinity by a smooth null hypersurface. If the 2-surface $ is a 
smooth spacelike cut of the future null cone of a point o e M, then the two independent point spinors at o 
can be transported to $ by the same law of transportation [15], too. Dougan and Mason [16] suggested to 
choose anti-holomorphic or holomorphic spinor fields, which were shown to form two complex dimensional 
vector spaces in the generic case. Their constructions are genuine quasi-local in the sense that they can be 
applied in any spacetime for generic spacelike 2-surfaces which are homeomorphic to topological 2-spheres. 
The properties of these suggestions have been studied in a number of situations [14-20]. In particular, 
the Dougan-Mason energy-momentum is null if and only if the domain of dependence -D(S) of a spacelike 
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hypersurface £ with smooth 2-boundary $ is a pp-wave geometry and the matter is pure radiation, and the 
energy-momentum is vanishing if and only if D(E) is flat [18-20]. Furthermore, they have been calculated for 
small spheres of radius r [15,17]. All the three expressions coincide in the leading order: in presence of matter 
they are ^nr 3 T ab tb, in complete agreement with the expectation, whilst in vacuum (or if at least an open 
neighbourhood of the point o is Ricci-flat) they are r 5 T abcd t b t c t d , where T a 

bed is the Bel-Robinson tensor 

and G is Newton's gravitational constant. The Dougan-Mason energy- momenta have also been calculated 
with accuracy r 6 , and the two constructions deviate in this order. 

Because of the lack of any geometric meaning of the coordinates, the notion of angular momentum in 
general relativity is a more delicate problem, and it is argued e.g. in [21] that angular momentum should 
be connected with the internal Lorcntz rotations of the theory, and hence should be analogous to the spin. 
In fact, general relativity has a Yang-Mills formulation (see e.g. [22,23]), in which Bramson derived the 
conserved Noether current corresponding to the Lorentz gauge symmetry and its superpotential 2-form [24] 
(and which has been rediscovered recently [25]). The Bramson spin superpotential 2-form, denoted here 
by w(X,n) a b, is well defined for any pair of spinor fields. Bramson used the asymptotic twistor equation 
to specify these spinor fields at future null infinity, and obtained an expression for the global spin-angular 
momentum by taking its integral on a spherical cut of null infinity. Thus to have a reasonable quasi-local 
spin-angular momentum associated with the 2-surface $, the two spinor fields and [ia must be spccifield 
there. For them it seems natural to use the spinor fields of the quasi-local energy-momentum. In fact, 
Ludvigsen and Vickers defined their quasi-local spin-angular momentum [14] as the integral of the Bramson 
superpotential with the spinor fields that they used in their energy-momentum expression. Recently it was 
suggested to use the anti-holomorphic or the holomorphic spinor fields in Bramson's superpotential, which 
yield genuine quasi-local spin-angular momentum expressions [20]. (Dougan and Mason suggested quasi- 
local angular momentum expressions, too, by using the spinor parts of the holmorphic or anti-holomorphic 
local twistors on $ in the Nester-Witten 2-form [16]. In the present paper, however, their suggestion will not 
be investigated.) However we think that it is not enough to give a definition, but it should be clarified in 
various situations whether the new angular momentum expression really has the expected properties of the 
angular momentum, or, more generally, whether Bramson's superpotential serves an appropriate framework 
for finding the (quasi-local) measure of the angular momentum of gravity. In fact, in [20] the energy- 
momentum and spin-angular momentum expressions based on the anti-holomorphic unprimed spinors were 
calculated for finite axi-symmetric pp-wave Cauchy developments, and the energy-momentum was shown 
to be an eigenvector of the spin-angular momentum tensor. Therefore the null energy-momentum and the 
Pauli-Lubanski spin are proportional, a reasonable property that is shared by zero-rest-mass radiative matter 
fields in Minkowski spacetimc. 

In the present paper the quasi-local Bramson spin-angular momenta are calculated for small spheres 
of radius r for the Ludvigsen- Vickers and for the holomorphic and anti-holomorphic spinor fields. Since 
the angular momentum is the momentum of the energy-momentum, the leading orders are expected to be 
higher than those for the energy-momentum, namely r 4 for non-vacuum and r 6 for vacuum spacetimes. In 
fact, for non- vacuum spacetimes any reasonable total angular momentum expression can be expected to yield 
(1.2) as the leading term, simply because the energy-momentum tensor plays the role of the source for the 
gravity. The spin-angular momentum is a part of the total angular momentum, thus the same r-dependence 
is expected for the leading orders of the spin-angular momentum too. The prototypes of the small sphere 
calculations are those of Horowitz and Schmidt [26] for the Hawking energy [27] and of Kelly, Tod and 
Woodhouse [28] for the Penrose mass [29]. Recently Brown, Lau and York [30] calculated the Brown- York 
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energy [31] for small spheres and increased the accuracy of some of the spin coefficients.* Although the order 
of accuracy of these calculations is enough to calculate the quasi-local energy-momentum even in order r 6 (for 
the Dougan-Mason energy- momentum see [17]), the order of accuracy must be increased for the calculation 
of angular momenta. All the spin coefficients are needed with accuracy r 3 . Thus first we review the main 
points of the framework of the calculations and present the spin coefficients and curvature components with 
r 3 accuracy. Section three is a quick review of the Nester-Witten and Bramson superpotentials, some of 
their properties and those forms of their integrals that we use. Since, as far as we know, the Ludvigsen- 
Vickers energy-momentum has not been calculated in vacuum in r 6 order, in Section 4 first we calculate 
this and compare with those of Dougan and Mason. Then the angular momenta will be calculated for the 
Ludvigsen-Vickers, the holomorphic and the anti-holomorphic spinors, both in non- vacuum and vacuum 
spacetimes. We will see that our expectation was, in fact, correct: In the non-vacuum case we recover (1.2) 
in all the three constructions, but with opposite sign for the anti-holomorphic spinors. In the vacuum case 
all the three construction give the same leading term. It is of order r 6 and the factor of proportionality is 
the contraction of the Bel Robinson tensor and an other average of the approximate Killing vector (1.1). 
The results will be discussed and summarized in Section 5. 

Our conventions and notations are mostly those of [32]. In particular, we use the abstract index nota- 
tions, and only the boldface and underlined indices take numerical values, e.g. a = 0, ...,3, A = 0, 1 and 
A = 0, 1. The spacetime signature is -2, the curvature- and Ricci tensors and the curvature scalar are defined 
by -R a bcd X b Y c Z d := Vy(V z A a ) - V z (VyA a ) - V [y , z] A a , R ab := R c acb and R := R ab g ab , respectively, 
and hence Einstein's equations take the form G a b '■= Rab — \Rg a b — ~8nGT ab . On the other hand we use 
the GHP formalism in its original form [33], and we refer e.g. to equation (2.21) of [33] as (GHP2.21). 

2. Small spheres 

First, mostly to fix the notations and ensure the coherence and readability of the paper, we summarize the 
geometric framework and philosophy of the small sphere calculations of [26] and especially of [28] . Since in 
this approximation the GHP equations have a hierarchical structure, they can be integrated with arbitrary 
high accuracy. In the second half of this section we integrate the GHP equations with the accuracy needed 
in Section 3, namely with r 3 . 

Let o € M be an arbitrary point, {t a , x a , y a , z a } be an orthonormal basis in T M with t a future pointing, 
and parametrize the future celestial sphere at o by l a {9,<t>) := t a + x a sin # cos </> + y a smdsm<fi + z a cosQ. 
They are precisely those null vectors at o whose scalar product with t° is one. Let Af a := dI + (o), the 'future 
null cone' of o, let l a denote the tangent of its null geodesic generators coinciding with l a (9,cj)) at o and 
let r be the affine parameter along the integral curves of l a , l a W a r — 1, and r(o) = 0. Then in an open 
neighbourhood of o the set Af a — {o} is a smooth null hypersurface and $ r := {p € J\f \r(p) — r} is a smooth 
spacelike 2-surface and homeomorphic to S 2 . $ r is called a small sphere of radius r. (r, 9, <f>), or equivalently 
(r, C, C), forms a coordinate system on the smooth part of Af , where £ := expi^cot |. Let us complete l a 
to a complex null tetrad {l a , n a , m a , m a } such that n a is orthogonal to, and the complex null vectors m a , 
m a are tangent to the spheres of constant radius % r . If e\ := {o A ,i A } is the corresponding normalized 
spinor dyad on Af a such that l a V a o B = 0, then, apart from constant transformations, the complex null 

* For small spheres the Hawking and the Brown- York energies give the same result in the leading orders, 
namely ^-r 3 T ab t a t b in non-vacuum and ■^Qr b T abc< it a t b t c t d in vacuum. Thus, assuming that the weak en- 
ergy condition is satisfied, the Ludvigsen-Vickers, the Dougan-Mason, the Hawking and the Brown- York 
energies are all positive for small spheres both in non-vacuum and vacuum spacetimes, which property can 
be interpreted as some form of the local positivity of the quasi-local energy. 
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tetrad (and hence the spinor dyad too) becomes fixed. The dual spin frame is := {— t a,oa}- The basis 
spinors {o A (r 1 C, (), i A (r, C, £)} (and the spinor components in this basis) have direction dependent limits at 
the vertex o. In fact, if £j£ := {0 A , I A } is the normalized spinor dyad at o associated with the orthonormal 
vector basis {t a 1 x a 1 y a , z a } (which we call Descartes spinors), then at o 



In the coordinate system (r, £, £) the covariant directional derivative operators take the form D := l a \7 a = 
(£p), S := m a X7 a = P(-^) + Q(-§^) an d 5 := m a V a , where P and Q are functions of the coordinates. The 
choice of the spinor dyad above yields the following relations on the GHP spin coefficients k = e = p — p = 
t + fit — fJ = p' — p' = t' — (3' + fj — 0, and expressions for the GHP differential operators: ~\Df = Df, 
df = (fi — (p — q)/3 — qT)f and d' f = (S + (p—q)f3—pf)f, where / has type (p,q). In flat spacetime the relevant 
spin coefficients and the functions P and Q are K = e = T = r' = a = a' = 0, p = — -, p' = 1 -, (3 = — L C, 
and P — -y=^(l + CO an d Q = 0. Thus in flat spacetime 5 reduces to ^qS, where S := (1 + CC)(J|0> 
and e.g. df reduces to i„0/, where df := 5f + ^(p - q)(f = ^(1 + C0 1_S ^((1 + CC) S /) with 
s := \(p-q), the spin weight of /. The Ricci identities (GHP2. 21-24) and the primed version of (GHP2.25) 
and (GHP2.26), and the commutators (GHP2.31) and (GHP2.32) show that in general the deviation of 
the non- vanishing spin coefficients from these values is of order 0(r), e.g. p — — £ + 0(r), a = 0(r), and 
similarly P = ^=^(1 + CO + 0(r) and Q = 0(r). Thus the flat space values of these quantities can be used 
as the initial data in finding the approximate solutions of accuracy r k , k > 1, of the GHP equations. The 
spin coefficients k! and e' and the operator \) do not play any role in the small sphere calculations; and, in 
fact, the geomtry of Af does not determine them. 

To find these solutions with accuracy r 3 , first let us expand the spin coefficients and the functions P and 
Q as their flat space value plus polinomials up to third order in r, e.g. p = —-+rp^ v> +r 2 p^ + r 3 p' 3 ^ +0(r 4 ) 
and P = y^(l + CO + rP' 1 ^ + r 2 P^ + r 3 P^ + <9(r 4 ), where the expansion coefficients are functions of 
C and C- Similarly, let us expand the curvature components retaining their first three non-trivial expansion 
coefficients, too, i.e. for example ip = ip^ + rtp^ + r 2 ip^ + 0(r 3 ). (At this point we should note that the 
different expansion coefficients have different (p,q) types. For example V>o°\ ipo^ an d are °f type (4,0), 
(5,1) and (6,2), respectively. Thus to save the extra care in using the GHP formalism we use the operators 
q5 and o^-) Then substituting the spin coefficients and curvature components into the GHP commutator 
(GHP2.31) (applied to the type (0,0) functions C a nd Q and the Ricci identities (GHP2. 22-24), we get 

P = ^(1 + CC) { \ + + ±r><® + r 3 ((^ f + ^V) + ±4®] } + 0(r% (2.2a) 

Q = ^ + CC){^4 0) + y/^ + ^(l^ + 4 2) )} + 0{r% (2.26) 

P = -\ + + \r 2 ^ + ^r 3 ((^) 2 + + 4°>$°>) + 0(A (2.3a) 

- = |r^ 0) + \r 2 ^ ] + ^r»(2$>4°> + ^) + 0{r% (2.36) 



T 



3 

+ 



V(2« +2(4 o } +^?Vr +9(*£> +4 2) )) +0(r% (2.3c) 



45 

Then from the commutator (GHP2.32) (applied either to C or to C) the spin coefficient [3 can be expressed by 
the expansion coefficients of the curvature and their S- and <5-derivatives. To express the derivatives of the 
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curvature components by the curvature components themselves, use the GHP Bianchi identity (GHP2.33). 
We get 



+ ^ 3 (^ 2) + ^(C^ - + -^M?> + + 3$V<°>)- (2.3d) 

- ^(C((^) a + - 2C^4 0) )) + 

Then the Ricci identity (GHP2.21) doesn't give anything new. To get the spin coefficients p' and a' let us 
use the primed version of (GHP2.25) and (GHP2.26). They will be expressions of the expansion coefficients 
of the curvature and of their S- and (^-derivatives. Some, but not all, of the derivatives of the curvature 
can be expressed by the curvature components themselves. First observe that the 5- and o^-derivatives 
of the zeroth order curvature components, e.g. of ip^ and of <f>Qi , can be calculated from their definition 
using (2.1) above. (Geometrically, these formulae are not Bianchi identities, because they don't contain the 
derivatives of the curvature itself.) As a result p'^ and ct'W become algebraic expressions of the curvature: 

P ' = Yr r + * + - ^ + 6A (0) ) + 0(r 2 ), (2.3e') 

a' = \r{^-^)+0{r"). (2.3/') 

To determine the higher order terms let us use the difference of the Bianchi identities (GHP2.34) and 
(GHP2.37), from which oSip^ and can be expressed by the curvature components and by o^n? an d 

o^io* j respectively. Then, using the definitions and (2.1) above, the remaining derivatives can be expressed 
by the O-components of the (first and second) derivatives of the Weyl spinor and by the spinor components 
of certain irreducible parts of the derivatives of the Ricci spinor. Since in general these formulae are rather 
complicated and we do not need the general expressions, we concentrate only on the vacuum case. In vacuum 
we get 



p / =^-^(4 0) +^? ) )-^(^ 1) +^ ) )- 

-r 3 (\(4 2) +$>) - ^» + l4°W?)+0(r% (2.3e») 

-'=^? ) +- 2 (^ ) -^ e (VeV5) ,) + 

+r3 (ro^ - ^(^/^o* - to^ ] + h^** - i ( ^ ?))2 ) + ° (r4) - (2 - 3/ " } 

Here we used the notation (V e r/0fc an d (V^Vni/Ok for the V e - and V( e Vj)-derivative of the Weyl spinor 
at o, respectively, contracted with k t A and (4 — k) o A spinors. The remaining Ricci identity, namely the 
primed version of (GHP2.21), and the Bianchi identity (GHP2.34) do not restrict these components of the 
derivatives of the Weyl spinor further. 

By calculating the determinant of the induced 2-metric on $ r , one gets the 2- area element on $ r . It is 
d$ r := -im [a m b] = -i(PP - QQ)' 1 ^ A d(. Then D(d$ r ) = -2pd$ r , and hence if d$ := lim r ^ (^d$ r ) = 
— 2i(l + CC) _2 dC A d(, the area element of the unit metric sphere, then 
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d$, = r\l- \r><$ \r*4® - ^ [4 0) 4> ] - <C) 2 + 9^o } ] + 0(r 5 ))d$. (2.4) 

Thus in non- vacuum spacetimes d$ r = r 2 d$ + 0(r 4 ), whilst in vacuum d$ r = r 2 d$ + 0(r e ). 

Finally, by t AA = ^o A 5 A + l a I a the primed spin vectors at o can be expressed by the unprimed ones: 
oa 1 = 2tA'ct- C an d t A — t A c oc- Thus, combining these with Lemma (4.15.86) of [33], we get 

jf OAl (C,O-"0AjC,O^;(C,C)^^ 

= 2 l t A , iCl ...t A[Cl t B '^ D \..t B '~ D ~ ^o Al ...o Ak o Dl ...o Dn L B K..L B ^L C K..L Ci m = (2.5) 

{ 0, otherwise. 

All the quasi-local energy-momentum and spin-angular momentum expressions for small spheres reduce to 
integrals of this type. 

3. The Nester— Witten and Bramson superpotentials 

First recall that for any pair of spinor fields X A and [i A the Nester- Witten and Bramson 2-forms are defined 

by 



u(X,fi) ab :=-(n A '\7 BB ,X A - ^b'Vaa'Ab), (3.1) 
w(X, n)ab ■= - iA (A Ms)£A'S' • (3.2) 

If $ is any orientable spacelike 2-surface in M then their integral on $ defines a Hermitian and a symmetric 
bilinear map, respectively, from C°°($, S A ), the (infinite dimensional) space of the smooth covariant spinor 
fields on $, to the field of complex numbers. Thus if A A , A = 0, 1, is any pair of spinor fields on $ then under 
the transformation of these spinor fields by constant SL(2, C) matrices the integrals 

PAB ' TG / U(AA ' XB )afc ' := 8^G / W(AA ' AB)ab (3 ' 3 ' 4) 

J $ J $ 

transform as a Hermitean and a symmetric spinor, respectively. Since the 2-forms (3.1) and (3.2) arc 
superpotentials for the energy-momentum and spin-angular momentum, for appropriately chosen spinor 
fields A A these integrals are intended to define the quasi-local energy-momentum and the (anti-self-dual) 
spin-angular momentum, associated with the 2-surface $, of the gravity plus matter system, respectively. 
Every two dimensional subspace of C°°($, S A ) with some SL(2, C) structure gives a potential definition for 
the quasi-local energy-momentum and spin-angular momentum. If one of the spinor fields in the arguments 
of the Nester-Witten 2-form, e.g. A^, is constant on $ in the sense m e V e X A = and fh e V e X A = (whenever 
the spacetime geometry is considerably restricted [18-20]), then u(X°, A A ) a t is exact, and hence the 00', 01' 
and 10' components of P AA are vanishing, and this P AA is null with respect to the SL(2, C) structure 
above. For this P AA and the spin-angular momentum tensor J AA BB := e AB J A B + e J we have 
Paa>J AA ' BB ' = (J 01 + J &1 ')P BB ' + {5*5?; J 00 + <5 B (5 B ' J ' ')^ 11 '- Thus the null P AA ' is an eigenvector of 
jAA bb JQ- joo j g van i sn i n g ; i. e . if an d only if the pull back of w(X°, X°) a b to $ is exact. However in general, 
without additional restrictions on A^ too, that is not exact. Mathematically, J AB is a measure on $ of the 
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non-intcgrability of the vector basis E AA , := A^A^;. If both A A are constant on $, e.g. the restriction to $ 
of the two covariantly constant spinor fields in Minkowski spacetime, then w(X , A B ) a {, are exact and hence 
J AB is vanishing. Thus J AA BB is a measure how much the actual vector basis E AA , is 'distorted' relative 
to the constant basis of the Minkowski spacetime. 

Next let us concentrate on the small spheres $ r of radius r about o, and denote the corresponding 
integrals (3.3) and (3.4) by P AB and J AB , respectively. Although these quantities can be calculated directly 
using the definitions (3.3) and (3.4), the formula (2.4) and the expansions of the spinor fields given explicitly 
in the next section, we can follow the philosophy of [26,15,17] by converting the 2-surface integrals into 
integrals on the light cone too: Since for the Ludvigsen-Vickers, the anti-holomorphic and the holomorphic 
spinor fields lim r ^ -P AB = and lim r ^ >/ AB = 0, by the Stokes theorem these integrals can be calculated 
as the integrals of the exterior derivative of the superpotential 2-forms on Af£, the portion of the null cone 
M between the vertex and $ r , too. To do the small sphere calculations in this way one must have a volume 
3-form on the null hypersurface Af . It is chosen to be e a bc '■= — i3!ri[ m6m c ], because, for the naturally 
defined real tri-vector e abc :— i3!/[°m b m c l on M — {o} and area 2-forms e a b '■= — i2m[ a fh b } on the spacelike 
2-surfaces $ r one has ^ k e a bc — — ^a&c> £ bc — l a £abc and e a bc — 3nr £(, c i, and hence the volume element 3-form 
on N is ^e a bc = \£[ab n c] = d$ r A dr. Then the integral (3.3) takes the form 



P r AB ' -^gJ^I^^'U = ^g[ l {(DXi)S(6X%)o B ' + {DX%)l B ' (5Xi)o A - 

-{DXji)o A {5X%)l B ' - {DX%)d B '{5X J X)i A + ( 3 - 5 ) 
+ (SX%)o A (8X%)o B ' - {5X%)o A {SX%)o B ' - ^X%G e f o F o F ,}d$ r/ dr' , 

where we used the Sparling equation. Similarly, (3.4) can be written as 

J ^ = ^rl v [a ^(A A ,A B ) bc] = 

° (3 6) 

= 8~7rG Jo j$ {^( A ^ A ^) ^ + iA ° B ) ^{^l)o A o B }d% rl dr' . 

For small enough r these integrals can be expanded as a power series of r. To calculate these integrals with 
accuracy r fe , by (2.4) we need to calculate their integrand with accuracy A k ~ 3 \ and hence, because of the 
operator D, we need to calculate the spinor fields with accuracy A k ~ 2 ^ . However, as we will see explicitly, 
to compute (3.5) with accuracy r k it will be enough to calculate the spinor fields with accuracy A k ^ . This 
is due to the special nature of the Nester-Witten form. 



4. Quasi-local spin-angular momenta for small spheres 
4.1 The Ludvigsen Vickers spinors 

First let us recall the definition of the Ludvigsen-Vickers spinors in the small sphere context [15]. They 
are the point spinors transported from o to $ r along N Q by the propagation laws (DXa)o a = DXo = 0, 
(SXa)o a = 8'Xq + pX\ = 0, where the spinor components are defined by X\oa — X la '■= A^e^ := Xa, 
and the two initial values for Ao at o are the 0-components of the Descartes spinors £ A = {— Ia,Oa}- 
Ag(0) := £ A o A (0) = i^2((l + CC)^ an d Aj(0) := £\o A (0) = i^2(l + CC)^- The first of the propagation 
laws implies that the spinor components A A are independent of the affine parameter r, thus A A on J\f are 
given explicitly by these expressions. (If a spinor has a name index, too, then its spinor components in 
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the spin frame {e^ } will be written as subscripts and its name index as a superscript. Thus, for example, 
A° is the 1-component of the zeroth spinor, A^ := X\l a , whilst Aq is the 0-component of the first spinor, 
Aq := \ A o A .) Since the Ludvigsen-Vickers spinors are completely determined by the initial values, i.e. 
the Descartes spinors £ A at o, the metric e AB on the space of point spinors at o determines an SL(2, C) 
structure on the space of the Ludvigsen-Vickers spinors on $ r , too. It is this SL(2, C) structure that should 
be used to define the quasi-local mass as the length of P AB . In general the pointwise scalar product of the 
Ludvigsen-Vickers spinors, e.g. that for the basis spinors e A A A B , is not constant on $ r . 

Following the general prescription of [26,28,15,17], let us expand the 1-components as a power series of 
r. But since the operator D in (3.5) and (3.6) reduces the power of r, to have r 6 accurate integrals, the 
spinor fields must be calculated with accuracy r 4 : A A = A A <°) +rA A W +r 2 A A < 2 > +r 3 A A ( 3 ) + r 4 A A < 4 ) +0(r 5 ). 
Substituting this expansion into the second of the propagation laws and using (2.2) and (2.3) we get 



A A(0) 


— a& A A — o<5A A — 


1 c\ A 


A A « 






A A(2) 


= ^°o ) A A(0) -^ 


ho A o + qVv a o > 


A A(3) 


= ^ooV (0) -^ 


, (1) A A + 1 7/5 (1) A A 

>10 A + 12^!' ' 


A A (4) 




) 2 +^r^) A " (0) - 




-\{^-\^- 





(4.1.1) 



10 g Vl' J J A n 



Thus A A (°) arc just the 1-components of the Descartes spinors: A A (°) = S k l a . Therefore the spinor 
components A A and A A (°) satisfy d\ A = 0, d\ A ^ + ^A A = and <9'A A (°) = 0, too. The pointwise 
symplectic scalar product of the basis Ludvigsen-Vickers spinors is e^A^A^ = e AB (l + \r 2 (^ + k r3< t>oo + 



3^ 4 [4o + tOAoi?) 2 + s^o Vo° } ] + ^(r 5 )); i.e. the natural S'L(2, C) structure cannot be realized by the 
pointwise scalar product even in vacuum spacetimes. 

Before calculating the Ludvigsen-Vickers spin-angular momentum, for the sake of completeness let 
us calculate their energy-momentum with accuracy r 6 in vacuum spacetimes. As a consequence of the 
propagation laws, all the terms except the fifth in the integrand on the right hand side of (3.5) vanish. But, 
because of od\ A — 0, the integrand will be of order r 2 , therefore to have r 6 accurate integral, by (2.4) it is 
enough to approximate d$ r by r 2 d$. Finally, using (2.5), we get 



P r AB ' = -^r 5 T a bcd t»t c t d £ A £% + -L r H*(V e T\ cd )t b t c t d £ A £% + 0(r 7 ), (4.1.2) 



where K AB ' := £ A £ B > € T*M may be interpreted as a translation at o in the AB' direction, and T abcd := 
ipABCD'ipA'B'C'D' , the Bel-Robinson tensor at o. Thus the Ludvigsen-Vickers energy-momentum coincides 
with the Dougan-Mason energy-momentum based on the holomorphic unprimed spinors [17] even in sixth 
order. It is known [17], on the other hand, that at future null infinity it is the expression based on the 
anti-holomorphic rather than the holomorphic unprimed spinors that coincides with the Bondi-Sachs (and 
hence the Ludvigsen-Vickers) energy-momentum. Therefore the Ludvigsen-Vickers energy-momentum is 
interpolating between the anti-holomorphic Dougan-Mason energy-momentum for large spheres near the 
future null infinity and the holomorphic Dougan-Mason energy-momentum for small spheres. 

Next let us calculate the Ludvigsen-Vickers spin-angular momentum in non-vacuum spacetimes with 
accuracy r 4 . Since by (4.1.1) the integrand on the right hand side of (3.6) is of order r, we may write 
d$ r = r 2 d$, and, using Einstein's equations, a straightforward calculation yields 
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JAB = 1 3 T a6 i" U AB d$ + ^5) = ^4^,^,^^^ <A £ B) + q^S^ (4 l3) 

where T ab is the energy-momentum tensor at the vertex o, and the symmetric complex matrix valued 1- 
form if AB is given by (1.1). Therefore the Ludvigsen-Vickers spin-angular momentum for small spheres in 
non-vacuum spacctimcs gives precisely the expected result (1.2): the leading order is r 4 , and the factor of 
proportionality is the contraction of the energy-momentum tensor and the average of the null tangent of the 
light cone of o and if AB , a 1-form field that can be interpreted as the approximate anti-self-dual rotation 
Killing 1-form that vanishes at o. We would like to stress that the energy-momentum tensor appears here in 
a rather non-trivial way, in contrast to the r 3 order calculations of the quasi-local energy-momentum, where 
the energy-momentum tensor was present explicitly in the exact expressions just because of the Sparling 
equation. The integrand of (3.6) is still of order r in vacuum, thus to calculate the Ludvigsen-Vickers spin- 
angular momentum in vacuum with accuracy r 6 it is still enough to write d$ r = r 2 d$ by (2.4). But then the 
integrand should be calculated with accuracy r 3 . A direct calculation yields 

J r AB ^TJT^ r5TahCd f Ubl c K£ B d$ + 0(r 7 ) = 

144ttG J $ ^ 414 ^ 

4 Jorp ,AA',BB',CC',D'E c DFc(A-cE) ± /iL7l 

= ^q T 1 AA'BB'CC'DD't t t t £ t E t p + <J(r ). 

Therefore in vacuum the leading term is of order r 6 , and the structure of this expression is the same that of 
(4.1.3) with the Bel-Robinson tensor of the gravitational 'field', instead of the energy- momentum tensor of 
the matter fields. 

4.2 The holomorphic spinors 

Recall that a spinor field Xa on a spacelike 2-surface $ is called holomorphic if fh h WbXA = 0, which in the 
GHP formalism is equivalent to d' X\ + cr'Xo = and d' Ao + pXi =0. If $ is homeomorphic to S 2 then there 
are at least two, and for metric spheres there are precisely two linearly independent holomorphic spinor fields 
A A , A = 0,1, for which e AB A A A B is always constant on $. Apart from certain exceptional 2-surfaces (e.g. 
future marginally trapped surfaces, i.e. for which the GHP spin coefficient p is zero) the two spinor fields 
can be chosen to be normalized to the Levi-Civita symbol: e AB X\X^ = e AB . For small perturbations of 
metric spheres, e.g. actually for small spheres, the number of the holomorphic spinor fields is still two, and 
they span the spin space at each point of $ (see e.g. [19]). Therefore the pointwise symplectic scalar product 
defines a natural SL(2, C) structure on the space of holomorphic spinor fields. 

Following the general prescription of the small sphere calculations let us expand the components of the 
holomorphic spinor fields: A^ = A a ^ + ••• + r 4 XA^ + 0(r 5 ). Then the condition of holomorphy above 
yields a hierarchical system of inhomogenious linear partial differential equations for these components: 



d'X ^-X^ = ^(c^X^+Dt-^+Ef 



(4.2.1) 



A A 

where C~,...,iJ~ are explicit expressions of the expansion coefficients of the GHP spin coefficients and the 
functions P and Q; and they are vanishing in flat spacetime. Its general solution can be written as the 
sum of two spinor fields. The first, as we will see explicitly, is determined completely by the zeroth order 
holomorphic spinor fields, which turns out to be a linear combination of the Descartes spinors at o, and, 
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in addition, its first order term is vanishing. Thus there are two such independent spinor fields, and they 
have the structure Xa = cia^a + t 2 Xa^ + r 3 A^ 3 ) + where \a^ 2 \ A^ 3 ), ■■■ are determined by a^E^ 
(see below). They are particular solutions of the inhomogenious system (4.2.1). The second is the general 
solution of the homogenious equations and turns out to be the sum of arbitrary complex combination of 
the Descartes spinors in each order, i.e. it has the form rA^S^ + ... + r 4 A^£^ + 0(r 5 ). Hence the 
r k accurate homogenious solutions form a 2k dimensional vector space. However, as noted in [17] (and 
realized first in the case of 2-surface twistors by Kelly, Tod and Woodhouse [28]), these 'spurious' solutions 
correspond to the lack of any canonical isomorphism between the space of holomorphic spinor fields on $ r 
and $ r < with different radii r and r'; or, in other words, between the space of holomorphic spinor fields on 
$ r and the space of point spinors at o. Hence they should be 'gauge solutions', and neither the integral 
of the Nester-Witten nor that of the Bramson 2-form should be sensitive to them. (This issue has not 
been discussed even for the Nester-Witten form.) To check this, let Xa be any spinor field whose zeroth 
and first order coefficients in its power series expansion are constant linear combinations of the Descartes 
spinors, A^ = cla^a + ^a^a + t 2 Xa^ + let ja — ^AaSai a ^th order gauge solution, k = 1,2, 
and determine the leading order in the integral of u(X,j) a b and of w(X,j) a b on $,.. (A^ has the structure of 
a particular solution for 6a = and of another gauge solution for «a = 0. However, we don't impose any 
further condition on Xa-) Then the leading orders in these integrals are 

/ _ f 0(r k+3 ) 

(To calculate these integrals explicitly we need to use the full (2.4).) Thus the expected leading terms, 
namely the r 3 and r 4 order terms in non-vacuum and the r 5 and r 6 terms in vacuum spacetimes for the 
integral of the Nester-Witten and Bramson 2-forms, respectively, are in fact not sensitive to the spurious 
solutions, i.e. they are really 'gauge solutions'. However, without further conditions on the spinor fields A ,4, 
the next orders are already sensitive to them. But since it is only the first order gauge solutions that yield 
ambiguity in the next non-trivial orders by (4.2.2), it seems natural to impose the gauge condition that the 
first order term in the spinor fields Xa be vanishing; i.e. we exclude first order gauge solutions by hand. 
Note that although in general, without additional structures, it is meaningless to speak about "how much 
homogenious solution is contained in the particular solution A^" , because of the specific structure of the 
particular solutions (given below) this gauge condition is well defined. With this gauge condition the first 
two non- vanishing orders of the energy-momentum and spin-angular momentum will be independent of the 
remaining gauge solutions. The gauge independence of higher order terms cannot be ensured by similar, 
additional gauge conditions. 

The quotient of the space of the holomorphic spinor fields on $ r and that of the gauge solutions form in 
fact a two complex dimensional vector space. Up to the second order in arbitrary spacetime we can choose 
the following particular solutions: 

A* = o A £% + r 2 (2« + *8V - (*g> + - ^°oVW + 0(r% 

i i i (4.2.3a) 

= ^ + r 2 ((^? + - ±«$V + i(2<> + i$?V )# + 0(r 3 ). 

In the quotient space these form a basis which is normalized in the sense e AB X^X^ = e AB (l + 0(r 3 )). In 
vacuum these reduce to Dougan's solutions [17], and substituting them into (3.5) we recover the Dougan- 
Mason energy-momentum with r 5 accuracy. (Having imposed the gauge condition above, by (4.2.3) one 



m non- vacuum, 
in vacuum: 

(4.2.2) 

in non- vacuum, 
in vacuum. 
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could compute the Dougan-Mason energy-momentum in non-vacuum with r 4 accuracy. However we don't 
expect anything interesting in this order.) Substituting (4.2.3a) into (3.6) for the self-dual spin-angular 
momentum we get the expected result (1.2). 

To calculate the holomorphic Dougan-Mason energy-momentum and the holomorphic spin-angular mo- 
mentum in r 6 order, we need to know the holomorphic spinor fields with r 3 and r 4 accuracy, respectively. 
In vacuum they are 

A- = o A £% - lr 4 + <W } + 3V4°W + + *l>? ) $ ) yei+ 

(4 2 36) 

Af = + Jr»$Vtf + ^r»#V^^ ^ 

+ ^ (-^W + + 3v4°W + + wpwy** + ° (r5) - 

These are normalized with accuracy r 4 : e AB X-^X^ = e AB (l + 0(r 5 )). Substituting its r 3 accurate part into 
(3.5) we recover the result of [17], i.e. (4.1.2). Finally, substituting (4.2.3b) into (3.6), we obtain (4.1.4). 
Therefore the Ludvigsen-Vickers and the holomorphic expressions yield the same results even in r 6 order 
both for the energy-momentum and the spin-angular momentum. 

4.3 The anti-holomorphic spinors 

A spinor field A^ on $ is called anti-holomorphic if m b Wb^A — 0, which in the GHP formalism is equivalent 
to d\\ + p'Xq = and <9Ao + a\\ = 0. The philosophy and the calculations are quite similar to those in the 
holomorphic case, thus we present only the results. Apart from the 'spurious' solutions, for the independent, 
normalized anti-holomorphic spinor fields with accuracy r 2 in an arbitrary spacetime we can choose 

A A = o A £% + r 2 (VfV - (V4 0) - A<°V)# + «(r 3 ), 
Xt = L A Ejt + r 2 ((# A<°V - (4* - \* { $ \^)o A Ya + 0(r 3 ) 
Substituting these into (3.6) for the anti-self-dual spin-angular momentum we get 

J r AB = -^r 4 T AAIBB ,t AA 't B ' E e BF £ ( E A £f + 0(r% (4.3.2) 

i.e. that for the Ludvigsen-Vickers and the holomorphic spinors but with opposite sign. A pair of normalized 
holomorphic spinor fields in vacuum with accuracy r 4 is 
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(4.3.1a) 



A A = o A £% + r 2 (tf<V 4 0) o A )£% + ir 3 ((4 1} + 2t e (V e ^)i)^ - + 2t e (V e ^) 2 )o' 4 )^ + 
+ + 2/V + 4tHf) (V (e V^) 2 - + l(4 0) ) 2 )o A £t+ 

+ r\±(l^ + 2lf + 4iV)(V (e V /) ^) 1 - |4°Vf + ^fVf)^ + 0(r 5 ), 
= ^ + r\^\ A - J#?V)tf + 

+ ir 3 ((4 1} + 2i e (V e ^) 2 )^ - + 2t e (V e ^) 3 )o^ + J$ V)#+ 
+ r 4 (-l(«' + 2iV + 4i e i') (V (e V /) ^) 3 + ^^(V^fl,,- 

- iV^ 0) - ^ o, ^ o) + 3^ 0) ^ ?) - i^y £ A+ 

+ r*(^(«' + 2/f + 4iV)(V (e V /) ^) 2 - j^Vf + ^(4 0) ) 2 + ^ 0) ^)^ + 0{r% 

(4.3.16) 

Substituting its r 3 accurate part into (3.5) we recover the result of [17], namely (4.1.2) with the numerical 
coefficient ^ instead of in the second term. (The r 3 accurate solution wasn't given in [17]].) Finally, 
substituting (4.3. l.b) into (3.6) we get (4.1.4). 

5. Discussion and conclusions 

The integral J AB of the Bramson superpotential with the independent holomorphic or anti-holomorphic 
spinor fields on orientable 2-surfaces in arbitrary spacetimes and with the Ludvigsen-Vickers spinors in 
certain asymptotically flat spacetimes are well defined quasi-local observables of the gravity plus matter 
system. Their usefulness and interpretation as the spin-angular momentum, however, depends on their 
properties in specific situations, e.g. in the case of small spheres. To determine the general structure of 
the small sphere expression of such an integral one can follow the argument of Horowitz and Schmidt [26] , 
applied originally to the Hawking energy: A simple dimension analysis shows that the number of derivatives 
of the metric in the coefficient of the r k term in the power series expansion of the quasi-local spin-angular 
momentum expressions is (k — 2). But the lowest order tensorial quantities as expansion coefficients appear 
for k — 4 in non-vacuum and for k — 6 in vacuum spacetimes, and, apart from numerical coefficients, these 
are the energy-momentum and the Bel-Robinson tensor, respectively, contracted with the unit vector t a and 
the generator of the physical quantity in question. The numerical coefficients may be different in the different 
constructions. 

The present calculations for the specific constructions confirm this 'universal' structure of the small 
sphere expansions and the interpretation of J BB as the gravitational angular momentum. Moreover, 
apart from the sign difference between the Ludvigsen-Vickers and holomorphic constructions on the one hand 
and the anti-holomorphic on the other in non- vacuum spacetimes, the results in the leading order are exactly 
the same in all the three constructions, and it is the expected result (1.2). The spin-angular momenta coincide 
in vacuum in the sixth order, and the factor of proportionality is just the Bel-Robinson tensor contracted 
with the unit vector t a and an average of the approximate boost-rotation Killing vector that vanishes at the 
origin. We stress that the approximate boost-rotation Killing vector appears at the end of the calculations, 
that is not put in the exact formulae by hand, in contrast to the translations in the small sphere calculations 
of the quasi-local energy-momentum. This result shows up the proper interpretation of the Bel-Robinson 
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tensor as a quantity analogous to the energy-momentum tensor of the matter fields: it gives the main 
contribution to the gravitational self-energy-momentum, as has already been pointed out [26,28,15,17,30], 
and angular momentum at the quasi-local level. It might be worth noting that our previous result [20] that 
the null energy-momentum P AA ' of an (axi-symmetric) domain of dependence -D(S) is an eigenvector of 
the angular momentum tensor J AA BB can easily be recovered in the small sphere approximation: If S r 
is any smooth hypersurface whose boundary is $ r , the dominant energy condition is satisfied and P AA 
is null, e.g. P AA = S A S A P} 1 , then D(E r ) is known to be a pp-w&ve geometry and the matter is pure 
radiation [18,19]. Then by o e -D(£ r ) the energy-momentum and the Bel-Robinson tensor have the form 
Taa>bb> = M 2 A° \ a>\° b \%, and T AA , BB , CCIDD , = |</f A° A°', X° B X%, \%\%, A°,A£, , respectively. Now it is 
easy to compute P AA and, using (4.1.3) and (4.1.4), the angular momentum tensor J AA BB explicitly. 
Finally, one can check that P AA is, in fact, an eigenvector of J AA BB . The quasi-local Pauli-Lubanski 
spin vector associated with the small sphere $ r , defined by Saa' |£aa'bb'ccdd'P BB J CC DD , is 
however vanishing in the expected leading oders, i.e. in r 7 and r 11 , respectively. 

Since the angular momentum is not expected to have any definite sign (unlike the energy), the minus 
sign in (4.3.2) doesn't mean that the expression based on the anti-holomorphic spinors would be 'non- 
physical'. The relative sign difference between the holomorphic and anti-holomorphic expressions shows 
only that, roughly speaking, the vector basis built from the holomorphic and anti-holomorphic spin frames 
are distorted relative to the constant basis by the curvature just in the opposite direction in r 2 order. 

Having imposed the gauge condition one can calculate the r 5 and r 7 order corrections. In order r 5 
we don't expect anything interesting. The r 7 accurate calculations for the spin-angular momentum would 
however be more interesting, because the holomorphic and anti-holomorphic constructions are expected to 
be different, but probably the Ludvigsen-Vickers and the holomorphic constructions are still coincide in this 
order. But these calculations would be much more complicated because the spin coefficients and the spinor 
fields would have to be calculated with accuracy r 4 and r 5 , respectively. 

In non-vacuum the result, apart from the sign in the anti-holomorphic case, is precisely the expected 
total quasi-local angular momentum of the matter fields. Therefore it is natural to think of the integral J AB 
with appropriately chosen spinor fields as the sum of the total angular momentum of the matter fields and 
the angular momentum of the gravity itself. But then the question arises whether that represents the total 
angular momentum of the gravitational field too, as claimed by some authors, or, in accordance with the 
expectation of Bramson, only its spin part; and in the latter case whether J should be completed by some 
quasi-local orbital angular momentum or not. Although the results of the present paper show that in the small 
sphere context the Bramson superpotential serves an appropriate framework for defining the gravitational 
angular momentum, in other situations it may yield wrong results. A general analysis shows that such an 
orbital term would have to be the integral of the momentum of the Nester-Witten superpotential [34] , but it 
is not quite clear how this moment would have to be defined. To answer these questions one should consider 
other situations, e.g. large spheres near spacelike and/or null infinity, stationary axisymmetric systems etc, 
the subject of a future paper. 
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